Normal Modes of Black Hole Accretion Disks 



Manuel Ortega-Rodriguez^ 

Department of Applied Physics and Gravity Probe B, 

Stanford University, Stanford, CA 94305-4090; 

^ ; and Escuela de Fisica & Centro de Investigaciones Geofisicas, 

O \ Universidad de Costa Rica, San Jose, Costa Rica 

O . 
<N . 

^ ■ Alexander S. Silbergleit 

<C ; Gravity Probe B, HEPL, Stanford University, Stanford, CA 94305-4085 

O ■ 

cn : and 
CN ■ 

>■ ■ Robert V. Wagoner^ 

^ I Department of Physics and 

^ \ Kavli Institute for Particle Astrophysics and Cosmology 

\0 \ Stanford University, Stanford, CA 94305-4060 

O ■ 

Ph; abstract 

o ■ 
^ : 

c/3 . This paper studies the hydrodynamical problem of normal modes of small adi- 

^ . abatic oscillations of relativistic barotropic thin accretion disks around black holes 

■ (and compact weakly magnetic neutron stars). Employing WKB techniques, we 

' obtain the eigenfrequencies and eigenfunctions of the modes for different values 

■ " " ' of the mass and angular momentum of the central black hole. We discuss the 

properties of the various types of modes and show for the first time that modes 
covering the whole extension of the disk (full-disk p-modes) can exist within the 
studied thin disk model. However, these modes have relatively short wavelengths. 

Subject headings: accretion, accretion disks — black hole physics — gravitation 
— hydrodynamics — relativity 



^mortega@cariari.ucr.ac.cr 
^ gleit @relgyro . Stanford . edu 
■^wagonerQstanford.edu 



- 2 - 



Introduction 



For the last decade our group has been studying the hydrodynamical problem of the 
perturbations of equilibrium models of relativistic accretion disks through a normal-mode 
fo rmalism. This approach , which we have called diskoseismology, was inspired by the work 
of iKato and Fukud (Il980l ). who showed that the effects of general relativity can trap modes 



near the inner edge of an accretion disk around black holes as well as weakly magnetized 
neutron stars that are sufficiently compact. The importance of these hydrodynamical oscil- 
lation modes lies in the fact that they provide a potentially powerful probe of both strong 
gravitational fields and the physics of accretion disks, since most of them do not exist in 
Newtonian gravity. Their frequencies depend upon the angular momentum as well as the 
mass of the black hole. (Although we shall not explicitly consider neutron stars here, the 
results obtained will also apply to them to first order in the dimensionless angular momen- 
tum parameter a = cJ/GM'^, since their exterior metric is identical to that of a black hole 
to that order. This assumes that the innermost stable orbit of the accretion disk lies outside 
the neutron star and its magnetosphere.) 

While other appro aches t o the sub ject involve more local (e .g., dispersion relation) 



considerations [see, e.g.. iKatd (1200 ll ) and 



Lubow and Ogilvid (1l998l ) and references therein] 



our focus is instead on global modes, in which all particles oscillate with the same frequency. 
All of the modes we study are adiabatic, with discrete real eigenfrequencies. However, a 
discussion of the effects of visc osity is presented at the end of th e pa per. For review s of 
relativistic diskoseismology, see iKato. Fukue. and Mineshigd (Il998l ) and IWagonerl (119991 ). 



The fully relativistic stu dy of normal mo des has revealed the existe nce of three basic 

types of oscillation: g-modes (jPerez et al.lll997l ) (hereafter RDl), c-modes (ISilbergleit. Wagoner, and Orteg 

200ll ) (hereafter RD2), and p-modes. The p-modes are pressure-driven oscillations in which 

(unlike the g-modes) the pressure restoring forces act with the gravitational rest oring forces, 



makin g the corotating eigenfrequencies larger than the radial epicyclic frequency. lOrtega-Rodrfguez. Silbere 
( 120021 ) (hereafter RD3) studied the properties of fundamental p-modes with low frequency 
(those axisymmetric modes with no nodes in their vertical distribution, and with the small- 
est possible eigenfrequency) . A summary of our previous results for all three types of modes 



and their potenti al relevance in determining the mass and angu 



was presented by IWagoner. Silbergleit. and Ortega- RodriguezI (120011 ) . 



ar momentum of black holes 



In the present paper we aim at preparing all the tools needed for the a more general 
hydrodynamical analysis by means of an extended WKB method. They can then be applied 
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to studying modes of all types; in this paper, the emphasis is on p-modes, since the WKB 
technique has not yet been fully used to study them. We investigate various qualitative and 
quantitative properties of the modes which are axisymmetric as well as non-axisymmetric. 
They are located near the inner edge of the disk (p- and c-modes), close to but not including 
the inner edge (g-modes), in the outer region of the disk (p-modes), as well as covering the 
whole disk (p-modes). 

In Section 2, the basic formalism that governs our approach is summarized, and the clas- 
sification of modes is given, including full-disk p-modes never considered before. In Section 
3 the WKB solutions of the radial eigenvalue problem are obtained for modes of all types and 
locations. Section 4 deals with certain general properties of the vertical eigenvalue problem; 
in particular, it is shown that the number of negative eigenvalues (needed for p-modes) is at 
most finite, in the proper parameter range (the complete proof is found in the Appendix). 
The WKB solutions of the vertical eigenvalue problem are given in Section 5. Based on them, 
the non-existence of both inner and outer high-frequency p-modes is demonstrated; on the 
other hand, explicit asymptotic solutions for the vertical eigenvalues of g-modes and full-disk 
p-modes are found. Section 6 contains two applications of the developed techniques, namely, 
explicit asymptotic formulas for the eigenfrequencies of g-modes and full-disk p-modes are 
derived. In Section 7, a carefuU des cript ion of the unperturbed disk model, based on the 
results of Novikov and Thorne Jl973h and lPage and Thornd (Il974j ). is given. The model can 



be used in various studies of accretion disks. Section 8 describes how numerical results are 
obtained from the formalism, before Section 9 closes the paper with a general discussion. 



2. Basic Assumptions and Equations 

2.1. Hydrodynamic structure of the unperturbed accretion disk 

We take c = 1, and express all distances in units of GMjc? and all frequencies in 
units of (? jGM (where M is the mass of the central body) unless otherwise indicated. The 
unperturbed disk model and equations for perturbations are exactly the same as in our 
previous papers RDl, RD2, and RD3, so the account here is very brief. 

We employ the Kerr metric to study a thin accretion disk, neglecting its self-gravity. 
The stationary [djdt = 0), symmetric about the midplane z = 0, and axially symmetric 
{d/d(p=0) equilibrium disk is taken to be described by the standard relativistic thin disk 
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model (iNovikov and Thornelll973l : iPage and Thorndll974l ). The velocity components = 
= 0, and the disk semi-thickness h{r) ~ Cg/^l ^ r, where Cs{r,z) is the speed of sound. 
The key frequencies, associated with free-particle orbits, are 



Q{r) (1 -4a/r3/2 + 3aVr') 



1/2 



1/2 



(2-1) 



the ro tational, vertical epicyclic, and radial epicyclic frequencies, respectively ( lAliev and Galtsov 
198ll ). The angular momentum parameter a = cJ/GM"^ is less than unity in absolute value. 



The inner edge of the disk is at approximately the radius of the last stable free-particle 
circular orbit r = rj(a), where the epicyclic frequency /t(rj) = 0. This radius rj(a) is a 
decreasing function of a, from 1) = 9 through rj(0) = 6 to rj(l) = 1. So all the relations 
we use are for r > rj, where K{r) > 0. Note, in particular, that 



Q{r) > Q±{r) > K{r) , a > ; ^±{r) > Q{r) > K{r) , a < 



[2-2) 



For a = (a non-rotating black hole), Q{r) = Q±{r) > K,{r). We will also employ the outer 
disk radius, r^. For stellar mass black holes in low-mass X-ray binaries, < 10^ Mq/M (see 
RD3 for more details). 

We here consider barotropic disks [p = p{p), vanishing buoyancy frequency]. (A gener- 
alization to small non-zero buoyancy may be found in RD2, Section 4.3, and RDl.) In this 
case hydrostatic equilibrium provides the vertical density and pressure profiles 



p = Po{r){l-yy , p = po{r){l~yy-'\ g = l/{T-l)>0, 



(2-3) 



where F > 1 is the adiabatic index (for brevity and convenience, we will use parameters g 
and F alternatively). One has F = 4/3 within any radiation pressure dominated region of 
the disk, and F = 5/3 within any gas pressure dominated region. The disk surfaces are at 
y = ±1, with y related to the vertical coordinate z by 



h(r) 



F - 1 
2F 



{2gy/'c,{r,0) 



:2-4) 



and the disk semi-thickness h{r) specified by equation (2.19) of RDl. More information on 
the unperturbed disk is given in Section 7. 
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2.2. Hydrodynamic equations for the disk perturbations 

To investig ate the eigenmodes of the d isk oscillations, we apply the general relativistic 



formalism that 



Ipser and Lindblonj (Il992l ) developed for perturbations of purely rotating 
perfect fluids. Neglecting the self-gravity of the disk is usually a very good approximation, 
examined in Section 7. One can then express the Eulerian perturbations of all physical 
quantities through a single function 6V = 5p/{p[3uj) which satisfles a second-order partial 
differential equation (PDE). Due to the stationary and axisymmetric background, the angular 
and time dependences are factored out as 5V = K(r, exp[i(m0 + at)], where a is the 
eigenfrequency, and the PDE becomes: 

1 (^Z^dV\ ^ d_ f^dV_\ ^Pj^y^^_ .2-5) 
dr yuj"^ — dr J dz \uj'^ dz J 

Here [3 = dt/dr = U^, the mode corotation frequency 

uj = uj{r, cr) = o" + mfi(r) , 

and g"^^ = 1/grr is a. Kerr metric component in Boyer-Lindquist coordinates. As in RDl (at 
the end of section 2.1), the assumptions h/r <^ ^ r/h and \m\ <^r/h have been made 

here to simplify the equation, since h <^ r. The terms neglected do not change the results 
appreciably, except near the corotation resonance. 

Note that the flrst of these inequalities is violated in a small neighborhood of the coro- 
tation radius (rc), where uj{rc) = 0. Since Q{r) is a monotonically decreasing function, a 
corotation resonance exists in the disk (spanning radii ri < r < Tq) if the mode frequency 
I cr I is within the range 

\m\Q{ro) < |(t| < \m\Q{ri) . (2-6) 

Since p-modes are trapped where cu^ > k^, the resonance lies where they rapidly vanish. 
For g-modes (cj^ < i^"^), those of interest (i. e., excluding large values of |m| and the radial 
mode number n) have Vc < ri [thus violating the upper limit in equation (12-61) ]. as can be 
seen from Figure 3 and the Tables in RDl, as well as equation (16-51) below. 

Another assumption of strong variation of modes in the radial direction (characteristic 
radial wavelength = \V/{dV/dr)\ <^ r) ensures the approximate WKB separability of 
variables in this PDE, V{r,z) = Vr{r)Vy{r,y). The 'vertical' part, Vy, of the functional 
amplitude V{r, z) varies slowly with r, as do all unperturbed properties of the disk [such 
as the key frequencies ( 12-ip ]. Multiplying equation ( 12-5^ by {P'^pVrVy)'^c1{r,0) and using 
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equations fl2-3p and fl2-4p with cj. = dp/ dp gives 



d 



o?Vr dr 



dVr 



(u;^ — K^) dr 



d ( dV„. 



'^g^^pVy dy dy 



+ 



,2\ ' 



where 



and 



(2-7) 

(2-8) 
(2-9) 



a(r) = py/g^/cs{r,0) . 

Since the left-hand side of this equation is a rapidly varying function of r and the right-hand 
side is a rapidly varying function of y, they are both equal to a slowly varying separation 
function S{r). Then with p^^dp/dy = ~2gy{l — y"^)^^, defining the alternative separation 
function "^{r) = ujI{1 — S) gives the resulting ordinary differential equations for the vertical 
{Vy) and radial {Vj.) eigenfunctions: 

d^V„ 



dy"^ 



dV., 9 
dy 



d^Vr 



dr'^ {io'^ — K^) 



d_ 

dr 



dVr 

dr 



Vy = Q 



K = . 



(2-10) 
(2-11) 



Note that, according to equation (12-lip . the radial wavelength is approximately expressed 



as 



(2-12) 



As in RD2, we assume that 



a(r) = 7(r) 



{r - riY 



< < 1/2 , z/ > 



(2-13) 



where 7(r) is some function bounded from above and away from zero, varying slowly with 
radius, and tending at infinity to a limit 7(00) = > 0. The singularity at r = is due to 
the fact that the pressure at the inner edge vanishes in a typical disk model [p oc (r — r^)'^] 
if the torque does, in which case /i = (F — l)k/2r. Typically, k = 2 and F = 5/3, giving 



p = 2/5. If even a small torque is app 



ied to the disk at the inner edge, which is likely due 



to magnetic stress (iHawley and Krolikll200ll ). the pressure and the speed of sound become 
nonzero at r = r,. This corresponds to a nonsingular a{r) with p = 0. 

It is also true that z/ < 1 in most disk models. If it is less than one-half, some limitation 
on the disk size (given 700, and vice versa) applies, which follows from (12-91) and the thin 
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disk requirement of the previous section, that is, h{r) ~ Cs/f2 ~ 1/aQ ^ r. Details can be 
found in RD3. 

Together with the appropriate homogeneous boundary conditions (discussed by RDl 
and RD2, and below in sections 3 and 4), equations fl2-10p and (12- lip generate the vertical 
and radial eigenvalue problems, respectively, which are entangled through a and entering 
both of them. The radial boundary conditions depend on the type of mode and its capture 
zone, as discussed in the following sections. The coefficient a, the vertical eigenfunction Vy 
and eigenvalue (separation function) vary slowly with radius (comparable to the radial 
variation of the properties of the equilibrium model), as do uj{r) and ti;*(r). 

Our goal is to determine the vertical eigenvalues \&(r, a) and the corresponding spectrum 
of eigenfrequencies a from the two eigenvalue problems for equations (12-101) and (12-1 ip . 
Along with the angular mode number m, we employ j and n for the vertical and radial 
mode numbers (closely related to the number of nodes in the corresponding eigenfunction), 
respectively, as in our previous works cited above. It is worthwhile to note that the eigenvalue 
^ depends on the parameters r, cr, a, and m through one combination only, namely uj^, which 
is clearly seen from equation (12-lOp . 



2.3. Classes of modes 

2.3.1. Modes in the Lindhlad resonance frequency range 

The classification of the oscillation modes of accretion disks was explicitly discussed in 
RD2 and RD3, so here we just review it briefly, and make one addition that was missing 
before. Consider flrst the eigenfrequency range 

a_ < o- < (T+ , (2-14) 

where 

o"_i_(m, a)= max [— mfi(r) + /t(r)] , cr_(m, a)= min [— mf2(r) — K(r)] , (2-15) 

ri<r<ro ri<r<ro 

which is characterized by the existence of the Lindblad resonances, i.e., the two roots of the 
equation u;^(r) — K^(r) = 0, denoted r± = r±{a, m, a), ri < r^ < r^ < ro- 

P-modes are, deflned as those for which /ujI — 1 < 0; with a in this range, the 
coefficient in front of Vj. in equation (12- lip is positive within ri < r < r_ and r+ < r < 
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for p-modes. So, in the case fl2-14p two types of p-modes are possible: inner p-modes, with 
capture (oscillation) zone Vi < r < r_, and outer p-modes trapped within r_|_ < r < r^. The 
c-modes are those with the difference 1 — ^/cu^ changing its sign; as shown in RD2, their 
capture zone is Vi < r < Vc < r^, with Tc closely related to their eigenfrequency, which is the 
Lense — Thirring frequency at r = Tc. Finally, the g-modes are defined as those for which 
\E'/c(j^ — 1 > 0, which implies that their capture zone is between the Lindblad resonances r_ 
and r_|.. This completes the mode classification within the frequency range fl2-14p . 

Note that, as seen from equation (12-151) . 

(j±(0,a) = ±Kmax{o,) = niax fi;(r) , cr±(— m, a) = — (Jip(m, a) , 

ri<r<ro 

SO one can limit oneself to m > 0. Although a± = ±K,max for m = 0, for m > 1 both 
cr+ and cr_ are negative, because Q{r) > K,{r). Moreover, |cr+| is very small, which is clear 
immediately from its definition (12-151) . Indeed, for m > the quantity —mQ{r) + K,{r) is 
negative and goes rapidly to zero at large radii; therefore a+{m,a) = —mQ{ro) + n{ro) and 
|cT+(m, a)| <^ 1. The dependence of K^ax on a can be seen in the first row of Table 1 and 
also in Figure 5a of RDl. The values of — cr_ for m = 1, 2 are given in the second and third 
rows of Table 1. (The values of |o"+| ~ 10^^ for m > 2 and they are even smaller for m = 1.) 
The corresponding values of the radii for the quantities in Table 1 are found in Table 2. 



2.3.2. Full-disk p-modes 
In the eigenfrequency range complement to (12-141) . namely, 

(7 > (j+ or cr < a_ , (2-16) 

one has 

[^^(r)! = \a + mQ{r)\ > K{r) , 

hence (u'^ — k^) > throughout the disk, and no g-modes are possible; whether c-modes exist 
in this range is not yet established, although their existnece seems rather improbable. Still, 
the whole disk can be a capture zone for some p-modes, which we call full-disk p-modes. 
These modes seem to skip the investigator attention before. 

It is natural to discuss also the full-disk modes whose eigenfrequencies lie within the 'g- 
mode' range (I2-I4p : they clearly do exist. This is seen immediately from the radial equation 
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(12-1 ip . or, even better, from its transformed equivalent given by equation (13-31) below, by 
way of analogy with quantum mechanics problems. 

However, the absolute majority of such full-disk p-modes have a very specific radial 
distribution. Namely, they are strongly concentrated at either the inner or the outer region 
of the disk, and practically absent in the remaining part of it, which means that the modes 
reduce, in fact, to either inner or outer modes (studied by RD3 in the low frequency limit). 
Evidently, the full-disk modes with eigenfrequencies in the range (12-141) oscillate in both 
the inner, < r < r_, and outer, r+ < r < r^, parts of the disk, and tunnel through its 
central region, r_ < r < r_|_. The length of the latter is huge, unless the mode frequency 
is very close to either 0"+ — or o'_ + 0. Therefore, if (for instance) the amplitudes of the 
growing and decreasing (through r_ < r < r+) components of some full-disk p-mode are of 
the same order at r = r_, then the total mode amplitude is enormously larger in the outer 
region compared to the rest of the disk, so that effectively this p-mode turns out to be an 
outer one. In the opposite case, when the growing and decreasing component amplitudes 
are comparable at r = r+, one gets effectively an inner p-mode. All other cases are in fact 
non-generic because they require fine tuning of parameters. 



3. WKB Solution of the Radial Eigenvalue Problem 

In this section we construct radial WKB solutions for p- and g-modes. The detailed 
derivation of the corresponding solution for c-modes was given in RD2, section 5. In fact, the 
radial WKB equation for g-modes was also obtained before, yet in RDl; however, there it 
was derived from an approximate radial equation, instead of using the exact transformation 
of section 3.1 first introduced in RD3. 

The essential assumption for the validity of the WKB solutions is that the corresponding 
capture zone does not shrink, so that we can well separate the vicinities of the boundary 
points from the core of the domain where the standard WKB approximation works. This 
assumption fails only for low-frequency p-modes with m = 0; low-frequency fundamental 
p-modes, for which, in addition, j = 0, have already been studied by RD3. 
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3.1. Transformed radial equation and boundary conditions 

We will exploit the radial equation (12- lip in a convenient form derived in RD3. Namely, 
we replace Vr with W according to 

W- ' Vr - ^ dW ^ ^3_^^ 

(c(j2 — K^) dr ^ (1 — ^ /ojI) dr ' 

where the second equality is nothing else as the radial equation (12-111) itself. We then 
differentiate it in r, and change the radial variable r to 



r 

2/ 



T = a {r ) 



dr , r^, = const (3-2) 



(the absolute value here serves conviniently both the p- and g-mode cases). In this way, we 
arrive at the following equation which is to be solved instead of the equation (12- lip : 

d^W ij)"^ - K-2 

— ^P(r)W = ,, PM.-,^^-^ (3-3) 

For the eigenfrequency in the range (12-141) . the coefficient P(r) is positive in the inner and 
outer p-mode capture zones, < r < r_ and r+ < r < Tq, whose boundaries r± are the 
turning points of equation (l3-3p . For g-modes the capture zone is, as said, r_ < r < r^., 
with turning points forming its both boundaries. If, on the other hand, the eigenfrequency 
satisfies inequalities (I2-I6p . then P{t) is positive throughout the whole disk rj < r < Tq, 
being the capture zone of the full-disk p-modes. Expressed via the r variable, these capture 
zones are Ti < t < t_, < t < To, t_ < t < t^, and Ti < t < To, respectively; we use the 
natural notation tx = T{rx) for A = ±, i, o. 

The proper boundary conditions at the disk edges were formulated in RD3 as 

dW 



uj (rj) cos 9iW + sin 9i 



dT 



[cj^(ro) - K^{ro)] cos 9o W + sin 



= , (3-4) 
dW 



dT 



(3-5) 



with — 7r/2 < 6i o < 7r/2 parametrizing our lack of knowledge about the physical conditions 
at the disk edges. 

Equations (13-41) . (13-51) form the complete set of boundary conditions for the radial equa- 
tion (l3-3p in the case of the full-disk modes. However, for the inner and outer p-modes only 



- 11 - 



one of these conditions is employed, since the other boundary is formed by one of the turning 
points, either r_ or r+. In both cases we require the solution to decay outside the turning 
point boundary of the trapping zone, that is, for r > r_ in the inner and for r < r+ in the 
outer case. Consequently, neither of the conditions fl3-4p . fl3-5p is used at all for g-modes, 
which are required instead to decay both to the left of r_ and to the right of r+. The decay 
requirement concludes the formulation of the radial eigenvalue problem for modes of all types 
(including c-modes, see RD2). 



3.2. Solution for the outer p— modes 



For the outer radial problem the convenient variable is given by equation fl3-2p with 



T 



«2(r') 



r+ 



v[/(r ] 



dr 







j-To 


\ ^(^ ) " 









dr 



(3-6) 



The radial equation fl3-3p on the interval < r < Tq with the boundary condition (13-51) and 
the requirement for Wir^ to decay when r is negative forms the eigenvalue problem whose 
WKB solution is obtained in a standard way. Namely, for all relevant values of r except the 
vicinity of r = we use the WKB expression satisfying condition (l3-5p : 



cos 



P^l'^{r')dr' - arctan {k^ cot Qo) 



(3-7) 



r=ro 



Recall that in general u!^{r) depends also on a, m and a [via VL = Q{r, a) and Q± = fi_L(r, a)]. 
Also a{r) and the vertical eigenvalue generally depend on the adiabatic index F. This 
explains the dependence of k^ on the parameters. Note that, unlike the standard WKB 
situations, here we cannot automatically neglect the derivative of the amplitude P"^/'*(r), 
because, according to (I2-I3p . (l3-2p and (l3-3p . it grows as some positive power of its argument 
when the value of the radius, and hence of r, is large. However, the WKB analysis shows 
that this term can be neglected if 0.25 P~^^'^ dP/dr ^ 1 at r = Tq, which, by the equalities 
referred to above, reduces to 

a/nir,) > /i(r,)/r„ . (3-8) 



Since fi(ro) 



~ r. 



-3/2 



this condition is valid except perhaps for extremely small values of a, 



for which it should be specially checked. 
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Near the turning point r = 0, as usual, we set 

W oc Ai{-p^l^T), p+ = dP{0)/dr > , 

where Ai{z) is the Airy function of the first kind. The standard matching of this and 
expression (13-71) provides the balance of phases in the form: 




P^/^(r)rfr = 7r(n+ 1/4) + arctan(A;o cot 6^0) , n = 0,1,2,... . 



(for negative values of n the r.h.s. is negative, so they are discarded). Returning here to the 
integration over the radius according to formula (13-61) and using equation (l3-3p . we arrive at 
the desired equation for the eigenfrequencies: 




a(r)v/(l - ^/oj2)(u;2 _ ^2)^^ ^ 7r(n + 1/4) + arctan (fc^ cot 9o) , n = 0,1,2,... 



(3-9) 

This equation contains, in fact, two unknowns, the eigenfrequency a and the vertical eigen- 
value the second relation involving them is delivered by the vertical eigenvalue problem, 
sometimes in the form of an explicit expression for 1 — "^/ool, as in RD3 for the outer 
low-frequency fundamental p-modes. 



3.3. Solution for the inner p— modes 

In this case of the inner modes it is convenient to define r by choosing r^, = in equation 



dMD: 



r 



nr 













dr 



Ti 









\ ^(^ ) " 


/ «2(r') 









dr' . (3-10) 



The radial eigenvalue problem consists of the equation (l3-3p on the interval < r < r_ 
with the boundary condition (13-41) and the requirement for W{t) to decay when r > r_. As 
compared to that for the outer p-modes above, it is solved in a slightly different manner 
resembling the radial WKB solution for c-modes constructed in Section 5.2 of RD2. 

Again, near the turning point r_, 

W oc At[p^J^{T_ - r)], p_ = dP{T_)/dT < , 

and the WKB approximation for the core of the capture zone matching this is 



W oc P-^'Ht) 



cos 



P^I\T')dT' -Tl/A 



(3-11) 
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However, we cannot make expression (13- lip satisfy the boundary condition (l3-4p since the 
amphtude P~^/^(r) diverges at r = = by the force of equahties (12-130 . (l3-2p and (13-31) . 

Therefore we need one more approximation, for the right vicinity of r = 0. According 
to the equahties last referred to, equation (13-31) in this region becomes, to lowest order: 

-1/(1-2/.) 



^ 2M/(l-2M)p^ = 



a;2(r,)(l-2/i)2^/(i-2M) 



(1 - ^K^) 



a. 



(3-12) 

Its solution satisfying the boundary condition (l3-4p is expressed in terms of Bessel functions 



as 



W oc 



1-2/i 



h = 



A = s 



1/2 fl- 2^" 



1-/X 



(3-13) 



2(1-/.) ' 

Here T e{z) is the Euler gamma-function, and Jb{z) is the Bessel function of the first kind. 
Note that < 6 < 1/2 for < /. < 1/2. 

After matching the asymptotics of solution (13-130 at r — > -|-oo with that of solution 
(13- lip at r — > -|-0 we arrive at the eigenfrequency equation which, upon returning to the 
radius as the integration variable according to equation (13-101) . looks similar to the result 
(I33D: 



a(r)v/(cu2 - /t2)(l - ^i/ul)dr = 7r(n + 1/2) - $i 



n = 0,1,2,.. 



with 



$i = arctan 



7ib 



TEil-b) fx 



TE{l + b) \2 



2b 



sin 6j ± uj (r,) cos 6',' 



(3-14) 



(3-15) 



For the special cases in which 6i = 7r/2, 0, 



±- 



TT fl-2fi 



(3-16) 



4 V 1 -/i 

respectively. Again, the above should be solved (for a and \E') together with another equation 
arising from the vertical eigenvalue problem. Also, if applied to the inner low-frequency 
fundamental p-modes studied in RD3, equation (13-140 gives a result which is close to the 
one obtained there, despite the fact that the WKB approximation is not expected to be 
valid. 
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3.4. Solution for the full— disk p— modes 

Finally, we deal with the frequency range defined in formulas fl2-16p . for which both 
points r_, r_|_ are absent inside the disk. Accordingly, the radial equation (13-31) has no turning 
points in the trapping region, that is, in the whole disk. The boundary conditions are (13-41) 
and (13-51) . so now both boundary parameters 6i and 6^, are involved in the problem. 

Near the inner edge, exactly as in the previous section, we have solution (I3-I3p . In 
the remaining part of the disk the WKB solution satisfying the outer boundary condition is 
valid, 



W oc p-i/^(r) 



cos 



/To 
P^/^{T')dT' - arctan(fco cot Bo) 



(3-17) 



[ko is defined in (l3-7p ]. The standard matching of equations (I3-13P and (I3-17P provides, upon 
returning to the radial variable, the following eigenfrequency equation for full-disk p-modes 
[recall that a > cr+(m, a) or a < cr_(m, a), see Section 2.3]: 

/ a{r) y(^I^2rr;^2y(f3^/^^^ = ^(^^ + 1/4) _ $ . + arctan(A;o cot Oo) , (3-18) 
Jn 

with $j as in equation (13-150 . 



3.5. Solution for the g— modes 



The capture zone of g-modes is r_ < r < r_|_, so it is convenient to set r* = r_ in the 
variable definition (l3-2p , 



pr 


■ vl>(r ) 








1 









dr 



r_ = 



\E'(r ) 

^lij") 



dr . (3-19) 



The radial equation (13-30 on the interval < r < r+ with the requirement for W{t) 
to decay outwards of the capture zone forms yet another eigenvalue problem whose WKB 
solution is constructed similar to the above ones. Namely, for all values of r except those in 
the vicinity of the boundaries we use the standard WKB representation: 



W oc p-^'Hr) 



cos 



P^/\T')dT' - $3 



(3-20) 



with some $„. Near the turning point r = 0, as usual, we set 



W oc Ai{-p^_!\), 



dP{0)/dT < . 
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The standard matching of the two above expressions, similar to the one carried out in section 
3.2, shows that (n , n" below are integers) 

$g = n{n' - 1/4) . (3-21) 

Near the turning point r = r_|_ we set, accordingly, 

W oc Ai[p^_/^{T+ - r)], p+ = dP{T+)/dT > ; 

matching this with the expression f l3-19p in the limit r ^ r+ — 0, we find the balance of 
phases to be represented by 

r P^'^{T)dT -^g = Ti{n - 1/4) , 
Jo 

or, using the above value of 

/"^pi/2(r)rfr = 7r(n + l/2) , 

where we have denoted n +n' = n + 1 (this step is pretty similar to the matching of section 
3.3). Returning here to the integration over the radius according to formula f l3-6p and using 
equation (13-31) . we arrive at the desired equation for the eigenfrequencies: 

a{r)^{m/ujl - 1)(k2 - cj2)dr = Tx{n + 1/2) , n = 0,1,2,..., (3-22) 

with the range of the radial mode numbers n specified by the fact that the l.h.s. of this 
equation is positive. This is nothing else as the equation (4.1) from RDI, with the only 
difference that the arctangent present in the latter is replaced here by 7r/2. The difference 
stems form the fact that in RDI the eigenvalue problem was considered not over the capture 
zone r_ < r < r+, but on the whole disk rj < r < ^ oo. However, for all the eigenfrequen- 
cies found in RDI the argument of the arctangent is very large, so that its value practically 
coincides with tx/2. 

All radial WKB solutions, including the one for c-modes from RD2, section 5, are now 
ready to be used with the proper solutions of the vertical eigenvalue problem. The latter, as 
menitoned, can be found by any method appropriate in this or that case; in particular, the 
vertical WKB solutions found below in section 5 can be used. 
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4. Properties of the Vertical Eigenvalue Problem and p— modes 



4.1. Boundary conditions for the vertical eigenvalue problem 



We return now to the vertical equation fl2-10p and discuss, first of all, boundary condi- 
tions for it. In view of the symmetry conditions 



dV„ 



dy 







or 



y=0 







(4-1) 



for, respectively, even and odd mode distribution through the disk, we can consider only 
half, < y < 1, of the whole interval. Near the singular boundary point y = 1 equation 
fl2-10p can be written as 



9 



1 



+ . . . 



dVy 

dy 



+ 



9< 



Vy = Q 



(4-2) 



where the dots stand for regular terms [recall g = l/(r — 1) b y equa t ion 



23])]. Hence, 



according to the analytical theory of second order ODEs [see e. g. lOlverl (119821 )]. the general 
solution to equation fl2-10p near y = 1 has the form: 



Vy = CMy) + C2(i - yy''v2{y), g^i. 

Here Ci_2 are arbitrary constants, and f 1,2(2/) are certain converging power series of (1 
with vi 2(1) = 1. 



(4-3) 

y) 



For the most physically relevant range of values l<r<2 {g > 1), the second term 
in the solution (I4-3P goes to infinity at the boundary. Therefore in this case the requirement 
that Vy{y) be finite at y = 1 implies C2 = and specifies thus a unique, up to a constant 
factor, particular solution of the vertical equation fl2-10p . Hence, it proves to be a proper 
boundary condition for the vertical eigenvalue problem. 

However, no fundamental principle forbids one from having F > 2 (0 < (7 < 1). In 
such case the term with V2{y) is also finite at the boundary y = 1 (and even turns to zero 
there). Therefore the requirement |K/(1)| < 00 does not select one of the two particular 
solutions in the representation fl4-3p . hence it cannot serve as a boundary condition to the 
vertical equation (12-lOp . But, since 1 — g < 1, the derivative of this term at ?/ = 1 is infinite, 
and we again have to get rid of it, because an infinite derivative is physically unacceptable 
(see also RDl, Section 2.3). 



Finally, for F = 2 [g = 1), which is the only remaining point of the physic a l rang e 
F > 1, the second particular solution has a logarithmic singularity at y = 1 Olverl (119821 )]. 
and thus once again should be discarded. 



The above argument shows that for any F > 1 the correct boundary condition for the 
vertical equation (12-101) is that both the eigenfunction and its derivative be finite at the 
boundary, or, equivalently, that Vy{y) is regular at y = ±1. Note that, as observed in RDl 
and implied by equation (14-21) and the boundary condition, the vertical eigenfunction near 
the boundary can be written in terms of the Bessel function of the first kind as 



Vy{y) cx (1 - y)-^J2^i2\uj^Wgil-y)) ^ = (g - l)/2 = (2 - F)/2 (F - 1) . (4-4) 



Note that the physical meaning of the above mathematical conditions is that the La- 
grangian pressure perturbation should vanish at the disk surface, and its ratio to the unper- 
turbed pressure should at least remain finite there. (These are standard boundary condi- 
tions.) The complete demonstration of this can be found in section 2.3 of RDl. 

4.2. Lower bound for ^/tu^: the number of eigenvalues corresponding 
to p-modes is at most finite for 1 < F < 2 

A new fundamental fact about the vertical eigenvalue problem is that, given any uj^ > 
and 1 < F < 2 {g > 1), the spectrum of eigenvalues \E' is bounded from below. Namely, in 
this parameter range there exists a bounded function A{g,ujl) such that for any eigenvalue 
the following inequality is valid: 



The quantity A is positive, so p-modes, of course, can exist. However, since the spectrum 
of eigenvalues is discrete, only a finite number of eigenvalues for which ^^/cu^ — 1 < can be 
present. Accordingly, for any given cj* and 1 < F < 2, the number of vertical eigenfunctions 
corresponding to p-modes is at most finite. 

The proof of the inequality (I4-5P and the definition of A are given in the Appendix. 
Since the exact expression for A is too complicated for quantitative estimates, a relaxed but 
explicit inequality is also derived there, namely. 




(4-5) 




g [co^^ + ig~ l)/2]' 
2{9 - lyool 



+ 1 , ^ > 1 (1 < F < 2) . 



(4-6) 
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For large values of the magnitude of the r.h.s. grows as uj^, giving thus space for more 
negative eigenvalues [the same is true for the precise bound A{g,ujl) as well]. 

For r > 2, no analogs of the inequality fl4-5p are known, so far. However, the WKB 
approximations of the solutions constructed below indicate that the number of vertical p- 
modes is always finite, independent of the value of F > 1. 



4.3. Behavior of \E'(lj^)/u;^ 



For a mode of any kind the value of u"^ changes somewhat, and often significantly, 
through the mode capture zone. It is important to find out how the ratio \E'(c<j^)/c<j^, and 
thus the vertical structure of the mode, changes through it. 



To check on this, we rewrite the vertical equation fl2-10p in the self-adjoint form: 



d_ 
dy 



+2g KV+^(i-y^)] {i-yy-Wy = o 



(4-7) 



We then multiply this by Vy and integrate over the interval (—1, 1). Using integration by 
parts and the fact that Vy{±l) is finite, we obtain 



dVy 

dy 



dy 



and thus 



kV + VP (1 - y^)] (1 - yT' V; dy>(}, 



(4- 



which inequality is, of course, trivial when \E' > 0. 

We now assume that the value of ul in (14- 7p . and thus of and Vy, is slightly changed, 
and apply a standard perturbation procedure to determine those changes. In the first order 
of perturbations we arrive then at 



d 
dul 



!\ kV + ^(i-y^)](i-i/T"'K.'^z/ 



where the inequality holds due to 



1 < , (4-9) 

ut j\ {I - y^f V^^ dy 

. It shows that, for any values of parameters and any 



mode, the ratio '^{ujD/ujI is a decreasing function of ul. 

Let us now find out how ujI itself evolves through the disk, which behavior proves to be 
far from simplistic, depending on all the parameters (a, m and a) involved in 
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Note, in particular, that the ratio Q{r) /^^±{r) grows with r for a < 0, is unity for a = 0, and 
decreases for < a < aj_ = 0.953. For very rapidly rotating central objects with a± < a < 1 
the ratio is not monotonic, having a single minimum not far from the inner disk edge [see 
formulas (12-11) . (12-21) and RD2]. (For a > a±, Q± achieves a maximum value within the disk.) 
The ratio in question is close to unity for all values of a if r is sufficiently large. 

So, uj^ grows monotonically through the disk for m = if a < aj^. For m > 0, the same 
is true in the following cases: 1) a > 0, a < 0; 2) a > 0, not too small, and < a < a±] 
3) a < 0, large enough in magnitude, and — 1 < a < a±. Otherwise, ujl{r) is not monotonic 
through the disk, which is especially clear when a is in the 'g-mode' range (12-141) . and the 
corotation resonance (where = u = 0) is present. 



5. WKB Solution of the Vertical Eigenvalue Problem 



We now prepare for the WKB solution of the vertical eigenvalue problem. Near the 
boundary y = 1 approximation (l4-4p will be used, with its large argument asymptotics 
given, to lowest order, by 



Vy{y) oc 



cos 



2\u,Wg{l-y)-{7i/2) {g-1/2) 



(1 - ?/)9/2-i/4 



(5-1) 



To construct approximations for the remaining part of the interval y G (0, 1), we recall the 
self-adjoint form (14-71) of the vertical equation, and carry out the same kind of change of 
variable that we used in the case of the radial equation to reduce it to the WKB form (I3-3P : 



i-y 



,2\ -9 



dy . 



(5-2) 



For g > 1 {1 < r < 2) the range of the new variable is the positive semi-axis, li < g < 
1 (r > 2), then < C < C with 



a/tT Fj; (1 - g) 
2 TE{3/2-g) 



(5-3) 



[We set C* = oo in the first case {g > 1), in order to have uniform notations.] Of course, the 
inverse function y = y{() is always well defined, and the vertical equation in terms of ( is 



dC 



+ Q{C)Vy = 0, Q{C)^2gool 



(1 - y^ 



2\2<?-l 



(5-4) 
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This equation calls for a WKB treatment if tu^ ^ 1, or |\E'| ^ 1, or both conditions 
hold simultaneously. However, the first of them fails at the Lindblad resonances r_ and r_|_ 
(when exist) since, by equations 02-11) and (12-21) . 

^*(r±) = «:2(r±)/fii(r±) < 1 . (5-5) 

This might not affect the validity of the WKB approximation for the relevant modes (inner 
and outer p-modes, g-modes), for instance, because the second condition might be still 
valid there, or for some other reason, suich as just the extended parameter validity typical 
for various WKB approximations. Note that, as far as the p-modes go, the second condition 
can apparently be valid only for those of them which have negative vertical eigenvalues. 



5.1. Solution for g— modes and p— modes with > 

5.1.1. Derivation of the WKB equation 

In this case Q is positive within the interval ?/ G [0, 1), or C G [0, C*), so, taking into 
account boundary conditions (14- ID . we can write the WKB approximation for the whole core 
of the disk as 

14ocg-^/^cos[0(?/)-T7r/2] , 



Jo 



i-y 



1/2 



,2\-l/2 



y 



) dy , (5-6) 



where T = 0, 1 for the even or odd eigenfunctions, respectively. For ?/ 1 — to lowest 
order we have 



Q^2''^gcol{l-yf3-\ and 0(y) ^ 0(1) - 2|cu,| ^^7(1 - y) 
so the corresponding asymptotic formula for the eigenfunction is 

[l)-2\uj,\y/g{l-y)-l7i/2 



cos 



Vy oc 



(1 - y)9/2-l/4 



(5-7) 



Matching this with the asymptotic expression (15-11) of the boundary solution (14-41) . we arrive 
at the desired equation for the positive eigenvalues: 



V ^~y^ ^ 



\uj:t\y/2g 



(5- 
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where j is any suitable integer, and 5r (introduced in RDl) is: 

5r = (3 - r)/4(r - 1) (even modes); 5r = (r + l)/4(r - 1) (odd modes) . (5-9) 



Equation (15-81) coincides with the equation (4.2) found in RDl for the zero buoyancy 
case, but now it is derived in a more accurate way, because an approximate WKB expression 
of the vertical eigenfunction for the core of the disk was exploited there. Instead, our formula 
(15-71) follows from the full equation (15-41) obtained by the substitution (15-21) . 

The l.h.s. of equation (15-81) can be expressed in terms of complete elliptic integrals. 



which proves to be very convenient 



integrals can be found in iDwight 



'or th e analysis; all the needed information on the elliptic 



fll96lh . 



5.1.2. Explicit asymptotic solution for g-modes 



By definition, \E'/co'^ > 1 for g-modes; due to this the WKB equation (15-81) can be esily 
written in terms of the complete elliptic integral of the second kind, E{k), as 



^E(v/l-cu2/^) = vr (j + 5r)/v^ • (5-10) 

This is exactly the equation (4.12) of RDl, which has been used there for the effective 
numerical solution. All the results thus obtained showed that, in fact, \E'/co'^ ^ 1 (see Tables 
1 — 3 in RDl). This allows for a simple explicit solution for \E': indeed, the last inequality 
implies the limit 

k=^/l- C^2/^ ^ 1 - , 



in which case, according to the expansion 774.3 from iDwightl (1l96ll ) 



E{k) = 1 + 0.5 In 4/A; -0.5 iky + 



k' = Vl - fc2 = \uj,\/V^ +0 



The substitution of this expression in equation (15-lOp provides the following asymptotic 
solution with the first order correction: 

■K {j + 6r) TT (j + 5r) 



2gE{k) 



29 



^9 

1-0.5 ( In 



1-0.5 



In [A/k 
- -0.5 



0.5 

'^9^1 



{k' 



(5-11) 



vr (j + 6r) ) 7r2 (j + bv)' 

This result without the correction in square brackets is given in RDl; with the correction, it 
reproduces numerical values of from Tables 1 — 3 in RDl very accurately. Note that the 
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above validity condition for the found solution converts, by the force of the solution itself, 
into the following inequality 

v^/l-*l - ^^1^ » 1 (5-12) 
v 2g\uo^\ 

Note that to lowest order ^ does not depend on the radius, i. e., on u;*(r). 



5.1.3. Solvability criterion and the non-existence of p-modes with > 

In the range < ^/cj^ < 1 relevant for the p-modes the WKB equation fl5-8p can be 
reduced to 

E(V^^^.)='Jl^, (5-13) 



\uJ^\^/2g 



with the help of the expressions 781.05 and 781.15 from iDwightl (Il96ll ) 



The minimum possible value, jmin, of the vertical mode number j is specified by the 
obvious requirement j + > 0. Since, for both mode parities, 6r decreases for F > 1, 
Jmm(r) is monotonically growing, i. e., larger values of F are more restrictive for the vertical 
mode number. The standard unperturbed disk model (see Section 9) employs F = 4/3 in 
the middle region of the disk and F = 5/3 in both inner and outer regions. Due to this, at 
least a part of the capture zone of any p-mode, be it an inner, outer, or full-disk one, lies 
in a domain with F = 5/3; thus for any p-mode of either parity, jmin = jmm(5/3) = [by 
expressions f l5-9l) . ^5/3 = 1/2, 1 for even and odd modes, respectively]. 

A solvability criterion for fl5-13p . 



'^,(r)| <j + 6r< ^Mr)\ , (5-14) 



follows from the inequality 1 = E{1) < E(k) < E{0) = 7t/2 valid for < < 1, since 
the l.h.s. of the equation decreases monotonically. Therefore, a p-mode with \E' > and 
j > jmin can exist if and only if inequalities fl5-14p hold throughout its whole capture zone. 
This means, in particular, that, for all the relevant radii, the number j + 6r must be within 
an interval of the r-dependent length (1/2 — 1/tt) ^/2g \uj^{r)\, which interval also moves 
along the r axis. It follows that equation f l5-13p has, in fact, no acceptable solutions at all 
in the whole disk. 



First, no inner and outer vertical p-modes with > are possible in the WKB approx- 
imation. This has been proved correct analytically [by examining the solvability criterion 



- 23 - 



f l5-14p at the Lindblad resonances and showing it fails at least at one of them], and then 
checking the result numerically. Seecond, the situation with the full-disk modes is not much 
better: one needs a very special combination of parameters, including the mode eigenfre- 
quency, for a single such mode to exist. In this case the reason is, essentially, that the mode's 
trapping zone - the whole accretion disk - is too large. 

To show the nonexistence, we note that no \E' > full-disk mode with eigenfrequency 
a and angular number m is present if, for at least one radius r < Tq, 

\u,{r)\<{2/'K)\uj,{ro)\. (5-15) 

Indeed, this inequality means that the intervals (15-141) taken at the outer edge and at some 
other point of the disk have an empty intersection. Clearly, (15-151) is true for m = 0, when 
it reduces to just ^i_{ro) ~ To < (2/7r)fi_|_(r). 

If m > 1, inequality (15-151) holds for both positive, a > > 0"+, and large enough 
negative, > cr_ > a, eigenfrequencies (see Section 2.3). Depending on the range of a, 
one takes either r = rj or r = r_|_ and exploits the fact that f2_|_(ro) is very small. Thus 
the only case when the full-disk > modes seem possible is the one with the negative 
and extremely small eigenfrequency, > cr ~ — mf2(ro) ~ —mro^^'^ > so that a;(ro) 
and ciJ*(ro) are both very close to zero. But then the applicability condition of the WKB 
approximation is violated. Therefore generically the full-disk p-modes with \E' > do not 
exist (within our WKB approximation). Our numerical results confirm this conclusion. 



5.2. Solution for p— modes with \E' < 

5.2.1. Derivation of the WKB equation 

In this case the coefficient Q of equation ( 15-4^ turns to zero at the point 

y+ = {l-co'J^)-'^' = {l + col/\^\)-'^' e (0, 1) , (5-16) 

with Q < for < ?/ < y+, and Q > for ?/_|_ < ?/ < 1; we write = C{y+), from equation 
(15-21) . The quantity y+ = y+{r, a, m, a) is a convenient parameter which we use freely in the 
sequel instead of and along with the eigenvalue "^/uj^ = — ?/^/(l — y^). 

Because of the turning point (+, we now need to explore four different approximations 
of the solution to the vertical equation (l5-4p . namely: 
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1) For < C < C+ and away from C+, the standard non-oscillating WKB approximation 
is used, 

VyOc\Q\-'/'coshE or |g|-i/SinhS, ^iC) = \QiC')\'^'dC' , (5-17) 

Jo 

for even and odd modes, respectively. When ( ^ — 0, its asymptotics is 

V, cc \Q\-'" {Aexp|-(2/3)Qf (C+ - C)'"l ± /l"' exp|(2/3)Qf (C+ - C)"']} ; (5-18) 



> 1 . (5-19) 



2) In the vicinity of the turning point C = C+ the standard approximation containing a 
linear combination of both Airy functions is used 

Vy cc Az[-QY\C - C+)] + KBz[-Qf{C - C+)] , 
with K =const. Its left far-field (C — C+ ~^ — oo) asymptotics , 



Vy oc (C+ - Cr'^' \ -exp[-(2/3)Qy^(C+ - O^^^'] + exp[(2/3)gf (C+ - C)=^/'] 



matches to fl5-18p giving 



K = ±1/(2A^ 



for the even and odd modes, respectively. The right far-field (C ~ C+ ~ 
formula then reads: 

Vy(x{C- C+)"^/^ COS [(2/3)gy^(C - C+f^^ - 37r/4 t arctan(2A2 



(5-20) 
-oo) asymptotic 

(5-21) 



3) For C+ < C < and away from both ends of this interval, the standard oscillating 



WKB approximation with some phase $ is used, Vy (x Q cos J^^ Q^^'^d(' + $ 
asymptotics at C C+ + is described by 



whose 



-1/4 



COS 



(2/3)gf (c-c+P + f 



Matching this with fl5-2ip provides the value of the phase: 

$ = -37r/4 =F arctan(2A^) . 



(5-22) 
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The behavior at the opposite end, C ^ C* ~ 0, is found to be 



Vy oc Q 



-1/4 



COS 



(5-23) 



4) Finally, near the right boundary C = C* (z/ = 1); solution (14-41) is used, as in Section 
6.1. Its left far-field asymptotics (15-11) . when matched with (I5-23P taking (l5-22p into account, 
leads to the desired eigenvalue equation: 



'l-(l-Vl/M2)(l-^2 



dy 



n 



(j + Ar 



y+ 



Ar = ig + l)/2 ± (l/vr) arctan(2A2) . 



(5-24) 
(5-25) 



Here the positive sign corresponds to even, and the negative sign to odd modes. The pa- 
rameter A defined in equation (I5-19P can be expressed in terms of complete elliptic integrals 
of the first and second kind: 



My+) = exp <; \uj^ 
[note that A ^ 1 — when 



29 {I -y"; 



-1/2 



[E{y, 



yl)K{y+)] 



(5-26) 



-0 (^ -0), and A +oo when ?/+ ^ 1-0 (^ 



-oo 



5.2.2. Solvability criterion and the non-existence of inner and outer p-modes with > 

As before, the l.h.s. of the eigenvalue equation (l5-24p can be conveniently written via 
elliptic integrals: 

E{,/l^)-ylK{y^l^) ^ vr (j + Ar) ^ _ 
a/1 - yl 1^*1 
where K{k) is the complete elliptic integral of the first kind. There is still a very important 
difference between this equation and equation (15-130 . Unlike the latter, here the l.h.s. can 
be infinitesimally small. Indeed, it decreases monotonically in the interval < y+ < 1 to 
zero at its right end, in partucular, 

k-^ [E{k) - (1 - e)K{k)] = {nk/A) [1 + e/8 + 0{k^)] , k ^ ; (5-28) 

the maximum value of the l.h.s. at y+ = is unity. 

Therefore, given a and m, the jth vertical p-mode with < exists if and only if, for 
all the pertinent radii, 

< J + Ar < \uJ*{r)\/7i . (5-29) 
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The value of jmin is apparently determined by the left inequality here. In the same way as for 
the modes with \1/ > 0, the strongest restriction comes from F = 5/3 and yields jmin = — 1 for 
even and jmin = for odd modes. Otherwise, requirement fl5-29l) is essentially less restrictive 
than fl5-14p for the modes with \1/ > 0. The reason is that all the intervals (15-291) . for any 
values of the radius, have a common left end at zero. Hence, as soon as u!^{r) ^ throughout 
the mode trapping zone [i.e., a ^ —mQ{r) there], the intersection of all the intervals (15-291) 
is not empty: it is the shortest interval corresponding to some radius r = r^m such that 
\i^*{rmin)\ is minimal. If this length satisfies 

\/2g \uj*{rmin)\/7^ > j + Ar , (5-30) 

then at least one vertical p-mode with \I' < (and j = jmin) exists. 

Nevertheless, it turns out there are no solutions corresponding to either inner or outer 
p-modes [governed by the WKB equation (15-271) ]. This is once again proved analytically 
same way as before, and then checked numerically. 



5.2.3. Explicit asymptotic solution for full-disk p-modes with \1/ < 

On the other hand, the existence of full-disk (WKB) p-modes with \I' < is immedi- 
ately seen. In fact, the existence criterion (15-291) requires only that |L(j^,(r)| be large enough 
throughout the capture zone, i.e., the whole accretion disk. But, since the eigenfrequencies 
in the full-disk mode range (12-161) are not bounded, this can be clearly achieved by choos- 
ing a sufficiently large value of a. This argument is supported by the following asymptotic 
solution of the vertical WKB equation (15-241) in the limit ujI^ 1 and —'^/uj'l ^ 1. 

We have shown in Section 4.3 that the ratio /ujI is a decreasing function of ujI. We 
will now explore further properties of this function based on the equation (15-241) . First, we 
define 

(j + Ar) 

and realize that 1 — '^/oo'l depends on F, j and uo^ = uj^:{a,m,a,r) only through the combi- 
nation cD*. (Even though, according to (15-251) . Ar depends formally on 1 — ^/tu^ through 
arctan [2yl^(?/+)], it turns out that in the range of interest A —>■ oo, arctan(2yl^) 7r/2, so 
no such dependence is left.) One can thus use the vertical WKB equation (15-241) to compute 
what we call a 'universal curve' (tD^ versus 1 — ^/cj^), which is shown in Figure 1. The curve 
does have the predicted monotonic behavior. 
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However, if both ul and {—"^/ujI) have large enoguh values, the vertical WKB equation 
fl5-24p can be solved asymptotically, to give the following two-term expression: 

1 - ^/ul = ul/K^^ + 1/4 = (7r^,/4)2 + 1/4, (5-32) 
^ (8/(7) (j + Ar)' = iS/g) [3 + {g + I ± l)/2f . (5-33) 

This two-term asymptotics is also plotted in Figure 1. One would expect the two curves, the 
exact universal curve and the asymptotic one, to be close for large enough values of 1 — \E'/u;^. 
However, they turn to be in a very good agreement within a remarkably wider range: even 
for the case —'^/uj'l = 1, (cD^. ^ 1.7), the two curves differ by less than 2 %. This means that 
we can use the analytical formula ( 15-32P practically for all the full-disk p-modes treated 
within the WKB approach, since it begins to fail 'much later' than the WKB approximation 
itself. 

5.2.4- Non-existence of p-modes with \E' changing sign 

One should mention here yet another possibility to find some p-modes, namely, those 
whose eigenvalue changes sign somewhere within the capture zone, so it satisfies equations 
f l5-13p or ( I5-27P on either side of the point r = r^, where it vanishes. (Note that the results of 
Section 4.3 show that such a change of sign can typically occur just once.) For inner or outer 
modes, the simultaneous validity of both solvability criteria at r = r^, restricts the solutions 
to the odd modes with j = — 1 only, and then the proper solvability criterion cannot be 
fulfilled at either r = r_ or r = r+. An analysis of the full-disk modes demonstrates that 
this is possible only for even modes with r = 5/3, j^l, and in a very narrow range of radii 
which may actually be in the F = 4/3 region. Thus at most a negligible number of p-modes 
with a changing sign of exist within the framework of our WKB approach. 

6. Applications: Explicit Formulas for the Eigenfrequencies 
of g-Modes and Full-Disk p-Modes with ^ < 

6.1. Eigenfrequencies of g— modes 

In RDI it was shown that, for given m > and the radial mode number n, the g-mode 
eigenfrequency amnj < tends to the frequency range boundary cr_ (m, a) = —cTm defined in 
fl2-15p when the vertical mode number j 00, i.e., when the eigenvalue \E'j becomes large. 
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In fact, all the eigenfequencies from RDI, Tables 1-3, have \E'/ci;^ ^ 1. In this limit, we can 
write the radial equation fl3-22p as 

a{r) {^! j / ujI - - uj'^)dr ^ a{r) \/'^ /uol \J - uj'^dr = Ti{n + 1/2) , 

and then use the explicit asymptotic solution for ^ found in section 5.1.2. To lowest order, 
formula fl5-lll) from there provides 

r = ^ /'!i±^y n = 0,l,...; j + 5r>0, (6-1) 

with the r.h.s. assumed small {n bounded). The assumption is exactly consistent with the 
studied limit, in which the capture zone is shrinking, with r_|_ and r_ both tending to merge 
at the point where mVL{r) + k{r) has its maximum, i.e., 

[mfi(r) + K{r)] = 0, mfi(rm) + ^(rm) = mfi(rm) + = -o"- {m, a) = ■ 

ar r=r,„ 

Because of this, to lowest order we can take the various quantities under the integral in the 
equation (16-11) to be: 

a(r) ^ a(r„) = a^, uj^{r) ^ a;*(r„) = /t(r„)/f2_L(r„) = K^/^im ; 
/t(r) - u{r) ^ /t(r„) - [amnj + mQ{rm)] ^ Km + (^m - mQ{rm)] = 2k„ ; (6-2) 
K{r) + uj{r) = K{r) + [amnj + mn{r)] ^ cr^ + amnj - A(r - r^)^ = A(r+ - r)(r - r_) , 

where 

1 rf2 

A = A(m, a) = -- — [mfi(r) + K{r)] ; (6-3) 

2 dr^ 



r=r„ 



note that the identity in (\Q-2h implies, in particular, that r+ — r_ = 2^ [am + o"mni ) / A. 

Using (l6-2p and (16-31) . with the same accuracy as above we find the following value of 
the integral (16-ip : 

r 41^;^* - *^v^ V (r,- r)(r- r,) dr = 

Jr_ ^^*(^) '^m Jr. 

^^y2fc:;:A(r+ - r^)^ .JT{y^)dt = (a^ + a„„,) . (6-4) 

Jo V 2KmA 

From formulas (16- ip and (l6-4p we obtain the desired explicit expression for the g-mode 
eigenfrequency amnj < 0: 



_ _ I , 2y/^^{n + 1/2) _ 2KmX n + 1/2 

~r O'mnj \^mnj\ ~ ^ / ■ , r \ '\ / ,t, ' \^ ^/ 

a„l2_LmVr(j + dr) V am^l±m 
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Up to a notation = 2fcmA, the last formula coincides with the result (4.8) of RDI, where, 
however, no explicit value of \E' was given. In addition, the current derivation provides clear 
asymptotic validity conditions of the result. 



6.2. Eigenfrequencies of full— disk p— modes 

6.2.1. Derivation of the eigenfrequency expressions 

We proceed here in the way similar to the one we followed in the previous case of g- 
modes. Assuming and — ^/cj^ large enough, we introduce the result fl5-32p to the proper 
radial WKB equation fl3-18p . which yields a relation: 



an 



Y(cj2 - Kynl){uj^/K^ + 1/A)dr = 7T{n + 1/4) - $i + arctan(A;„ cot 60) . (6-6) 



Since ojI is large, the term arctan(/co cot ^^o) 7r/2 by formula (13-71) . unless 60 = vr/2, in 
which case it vanishes; in any case, the term is within the range [— 7r/2, 7r/2]. The same is 
true for the phase shift $j, defined by the formulas (I3-I5p . (13-130 and (I3-I2p . which depends 
on the inner edge boundary condition parameter 6i. Hence these terms, as well as the 
additional 7r/4, are very small compared to vrn ^ 1 [from equation (l8-2p ]. For this reason, 
we neglect them in the following results, which become thus independent of the boundary 
conditions within the whole range of parameters 6i, Oq. 

Now, upon expanding the square root on the right of equation (I6-6P in inverse powers 
of ujI, the l.h.s. of equation (l6-6p becomes 



aVL 



dr . 



Therefore, neglecting the higher order terms, we obtain finally the following analytical ex- 
pressions for the eigenfrequencies: 

c^oni = \/ {rrn - hj) /hj ] (6-7) 
m |±y^e^ + {irn - hj - rn^hj) /m?Iij - e^j , (6-8) 



where 
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'3j 



— - rfr , J4,- = / ^ • (6-10) 



Of course, only those values of the mode numbers m, n, and j should be considered 
which give real values of eigenfrequencies in the proper frequency range (12-161) . 

It turns out, though, that all the mentioned requirements are not at all restrictive for the 
modes with the lowest possible radial numbers; the real restrictions stemming from the 
validity of our approach are discussed in Section 8. 

Note that the relative magnitudes of the integrals lij involved in equations ( 16-70 and 
(16-81) are governed mostly by the quantities fi±(r) and fi(r), which rapidly decrease away 
from the inner disk edge. One thus expects to be by far the largest of them, with 
coming second one or several orders of magnitude smaller, and the remaining two far behind 
the second one; this is perfectly confirmed by their values given in Table 6. Moreover, using 
the mean value theorem for integrals, one finds that there exists a radius rj such that 

e, = h,/h, = fi(r,) ^ rf'/' « 1, r, » r, > 1 . (6-11) 



The fact that the value rj is large enough follows from the radial behavior of the integrands; 

Id 

< 10~^, for To = 10~^, 100, and 20, respectively. 



so, the number ej should be small, and indeed Table 6 gives it at ~ 10 ^, ~ 10 ^, and 



6.2.2. Eigenfrequency dependence on the mode numbers m, j, n 

The above insights become particularly helpful for the discussion of the behavior of 
I'^mnjl given by formulas (16-71) . (16-81) . as functions of the mode numbers n, m and j, assuming 
that all the unperturbed disk parameters, starting with a, are fixed by some thin disk model, 
such as the one described in the next section. 

The first relevant property is rather straightforward: apparently, all eigenfrequencies 
grow when n increases (and the two other mode numbers stay fixed in the allowed range, as 
assumed in all other cases as well), as can be seen in Tables 7 and 8. The dependence on 
m is more tricky, although it is not very difficult to see that {c^njl always decreases when 
m(> 0) increases. As for |cr~„j|, it evidently grows with m if e| — hj/hj > 0. It also grows 
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in the opposite case when 

but starts decreasing when m? becomes larger than the number on the right of this inequahty. 

The dependence on the vertical mode number seems to be the most complicated until 
one notices that ej = hj/hj would be entirely independent of j if the parameter g = l/(r — 1) 
were constant along the disk radius. If this is the case, would also be constant, would 
be pulled outside the integrals, and would thus vanish from their ratio; see formulas ( l6-9p . 
Hence the dependence of ej on j is very weak at most, the same is true for the ratio I^^j/ Iij, 
and the trend is defined mostly by the ratio I^j/Iij, which appears to be roughly quadratic in 
j. The data of Tables 7 and 8 do not, however, confirm or disprove this observation, because 
the eigenfrequencies with different values of j in them have also very different values of n. 

Finally, one can notice that, by the definition fl5-33p . 

(even mode) = Kj (odd mode) . 

Since the approximations (16-71) . (16-81) . and, in fact, all the eigenfrequencies found from the 
radial and vertical WKB equations, depend on j through the combination Kj only, the same 
equality is true for them as well: 

c^mni-i (even mode) = (odd mode) , (6-12) 

provided both sides of it exist for a given set m, n, j. Thus within the WKB approximation, 
two modes can be excited at each eigenfrequency; one of them is even, the other is odd. 



7. Unperturbed Hydrodynamic Disk Model 



To carry out the hydrodynamical analysis of the accretion disk oscillations, one needs, 
in the first place, to formulate the assumptions about the unperturbed state. In addition 
to what has already been said in Section 2, we make all the relevant statements on the 
unperturbed disk model in this section. The results discussed in it were obtained by assuming 



that the opacity is mainly d ue to (optically thick) e 



for r < 4 X 10^{L/ LEddf^^ fishapiro and Teukolskv 



ectro n scattering. This condition holds 



19831 ). so our results at larger radii are 



only approximate. {LEdd oc M is the Eddington luminosity.) 
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It is notable that the interior structure of the (zero buoyancy) accretion disk enters our 
hydrodynamical formulation only through the index F and the function a{r), defined by 
equation fl2-9p . (This is not to be confused with the usual viscosity parameter, here denoted 
by a*). The relativistic factor 

f3^ = (1 + ar-^/^) [(1 - 3r^^ + 2ar-^/^) (1 - 2r-^ + aV-2)]-i/2 (y.^) 

that appears in the definition of a{r) is of order unity except for values of a close to 1 and 
simultaneously r close to ^ 1. The value of the coefficient f l7-ip at r = ri{a) as function 
of a is shown in Figure 2. 

The values for the sp eed of sound are obtained from the fully relativistic results of 
Novikov and Thorne 

Cs/c = 5.0 X lO-3(L/L^,drf)i/5(«,M/M0)-i/i%-i/5(a)r-9/2Og-^/5(a,r)r'-i/i°(a,r) 



[Q{a, r) + 5Q(a, r)Y^^ (gas pressure dominance, F = 5/3) , 
[Q(a, r) + 5Q{a, r)] (radiation pressure dominance, F = 4/3) 



where 



5Q(a, r) = {l- F,)B{ri, a)C'/\ri, a)T{r, a)B-\r, a)(l - 2r-' + 



a r 



-2U/2 



(7-2) 
(7-3) 

(7-4) 



In these expressions, t] (shown in Figure 3) is the 'efficiency factor' which relates the mass 
accretion rate to the luminosity: L = rjMc^. (It is the energy lost by a unit of rest-mass 
energy as it accretes from far away to the inner edge of the disk.) The constant F^: is the 
fraction of the inflowing speciflc angular momentum at that is absorbed by the black hole 
(so there is a torque at u nless F» = 1). The fu i iction s A, B, V, and S are relativistic 
correction factors deflned in 



Novikov and T 



r, while Q [deflned in iNovikov and Thornd (119731 ): 



lorne 



19731). which approach unity at large 



Page and Thornd (119741 )] is the factor 



proportional to the stress, which is responsible for making a{r) singular at rj when there is 
no torque present at the inner edge. In this paper we consider the case in which there is 
such a torque, and thus 6Q has been added to Q [see IPerezI (119931 )]. Expressions (I7-Ip - (I7-4I) 
completely specify the key function a(r), deflned by equation (l2-9p . in the domains of both 
gas and radiation pressure dominance discussed next. In particular, the two critical powers 
involved in it and used in calculations are fi = and u = 9/20. 

Since we deal with all kinds of r aodes in this paper, we n eed to consider the dependence of 
the adiabatic index F on the radius (IOrtega-Rodriguezll2002l ) . The regio ns where gas or radia- 



tion pressure dominates were obtained by making use of the solutions of 



Novikov and Thorne 
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( Il973l ). We define the transition radius as the radius at which the gas pressure equals the 
radiation pressure. The results of these considerations are shown in Table 3. Gas pressure 
dominates for ri < r < Ri and r > R2, while radiation pressure dominates for Ri < r < R2. 
In our calculations we use F = 5/3 and F = 4/3 for these two respective regions. 

Typical values for the disk density and central object mass make the neglect of self- 
gravity a common assumption. An order-of-magnitude analysis ( jOrtega-RodriguezI l2002l ) 
shows that self-gravitational effects can be ignored, at both the unperturbed or the perturbed 
level (Cowling approximation), whenever (now in normal units) 

« 7T77TT^ ■ (7-5) 



GM/c2 (M/M0)O-96 

Even though this estimate assumes a* = 0.1 and L = LEdd, condition (I7-5P is only slightly 
less restrictive for lower values of the luminosity. Thus, for stellar mass black holes the 
condition is satisfied by many orders of magnitude, while one needs to be more careful with 
supermassive black holes when studying the outer modes. A black hole mass of IO^Mq, for 
example, implies the constraint r <^ 2 x lO^GM/c^. 



8. Numerical Results for Full— Disk P— Modes 

At this point, we calculate some eigenfrequencies, aiming at those with the lowest mag- 
nitudes consistent with the appropriate restrictions discussed below. We compute them for 
different values of a, j, m, Tq, considering only the even modes [see relation (16-121) ]. We also 
consider full-disk p-modes with negative essentially the only p-modes to be captured by 
our approach. 

We make use of the unperturbed disk model described in the previous section, and choose 
the following ('typical') values for the parameters: L/LEdd = 0.1, a^^ = 0.1, M = lOM©, 
and F^, = 0.95. We also choose = 10^, typical of stellar mass black holes in low-mass 
X-ray binaries (LMXBs). In addition, we perform calculations corresponding to values of 
To = 100 and 20. Recall also that, as shown in Section 7.2, our results are robust regarding 
the boundary conditions at the disk edges. 

Using the fact that 1 — \I'/a;^ ~ 1^* [equation (l5-32p and Figure 1] and that \a\/il_i_{r) 1 
over most of the disk, one obtains the estimate 



r/Xr{r) ~ 10 



(8-1) 
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at radii r ^ Vi from equation (12-121) . using also equation fl7-2p . Since the integrand in 
equation fl3-18p is 1/Xr{r), the integral is dominated by its value near Tq. 

Note that the fundamental requirement allowing for our WKB separation of variables, 
r/\r{r) ^ 1, is thus enforced by its value near rj. From equations fl3-18p and fl8-ip . we then 
obtain the radial WKB requirement 

n > (r„/r,)3 . (8-2) 

This lower limit is seen to agree approximately with the minimum values of n found in Tables 
7 and 8, obtained from another WKB requirement discussed below. 

We turn now to the issue of the restrictions on the eigenfrequencies to be found. The 
first of them is 

1 - <il/ujl > 1 (8-3) 

(WKB parameter 'not small'), and our computations prove that this is the key limitation: 
all other conditions are satisfied if this one holds. 

Indeed, the restriction (18-30 . in view of relation (15-320 . implies that 

\uj,\ > y^OJbiq . (8-4) 

The relevant values of Kj for several lowest values of j are given in Table 5. The mode 
existence criterion (15-291) . 

\uj,\ >7r(j + l + (7/2)/v^, 

turns out to be valid whenever condition (l8-4p is true, as one concludes from Table 4. 

We begin the computations by calculating Iij, l2j, Isj, I^j, and ej for different values of 
j and a (Table 6). The value e = 2.04 x 10~® is the same for all the Tq = 10^ entries in Table 
6; for r„ = 100, it is 0.00171 < e < 0.00181, while for = 20, one finds 0.0283 < e < 0.0325. 

Tables 7 and 8 list the main results. In them, a is calculated taking into account the 
dominant condition (18-41) . which determines the minimum radial mode number n (and thus 
the lowest possible \a\) according to the formulas (I6-7P or (l6-8p . Namely, these lowest values 
are achieved when condition (l8-4p is valid all the way through the disk except just a single 
point r = Tc where the inequality (l8-4p turns into an equality; lower values of n (and |cr|) 
would give an interval of radii where (l8-4p is violated. The critical radius Vc is also listed in 
the tables; usually it coincides with the inner radius rj, or the transition radius Ri, at which 
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point the pressure becomes radiation dominated (see Section 8), and has a jump. In 
a few cases, the critical radius hes within the radiation-pressure dominated region; in such 
instances, Tc is indicated in the table by the appropriate factor multiplying Ri. 

It is important to mention that, as stated above, we might have missed eigenfrequencies 
with moderate magnitudes (corresponding to \u!^\ ~ 1) in our calculations. Those are the 
eigenfrequencies whose valus are too high to be captured by the low frequency analysis of 
RD3, and at the same time too low to be properly described by the current WKB approach. 
Those missing should include some inner and outer modes, as well as full-disk modes with 
radial mode numbers somewhat smaller than the ones found here. 

For instance, for m = and j = —1 (and gas pressure dominance) the WKB criterion 
(18-41) gives \a\ > {3/2)Q±. Now for most values of a, the maximum value of fi_L(r) > Sumax- 
As also seen from Tables 1 and 7, this WKB criterion then gives |cr| > 5Kmax, which should 
allow lower eigenfreque ncies that are greater than Km n ^ to e xist. The (m = 0) traveling 



acoustic waves found by lHonma. Matsumoto. and Katd (1l992l ) and lMilsom and TaamI (1l996 



19971 ) had a narrow range of frequencies near Kr 



9. Discussion 

As can be appreciated from the tables, for the hydrodynamical modes that we found 
the usual trend is to have an eigenfrequency which increases monotonically with a, except 
for a possible maximum for high values of a (like the low-frequency fundamental p-mode). 
The dimensional frequency / = 3.23 x 1O^(M0/M)|(t| Hz. For the mass of the central object 
that we adopted (M = IOMq), the minimum eigenfrequencies in Tables 7 and 8 are thus (in 
order) 378, 552, 568, 769, and 988 Hz; all but one corresponding to a = 0. Moreover, the 
eigenfrequencies increase monotonically with j. 

Comparing results for different values of the outer radius Vo, we find that the smaller the 
disk, the smaller the value of n for a given eigenfrequency. Also note that the value of the 
lowest eigenfrequency is independent of Tq. This is because the eigenfrequency is calculated 
using condition f l8-4p . 

It has been found that the pai r of high frequency QPOs seen in a few black hole LMXBs 



have a frequency ratio close to 3/2 (lAbramowicz and Kluzniak 



2001 



McClintock and Remillard 



20061 ). We note that the ratio of the j = 1 to the j = eigenfrequencies of the m = modes 
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is 1.57, independent of the value of a, as seen in Table 7 (for Tq = 10^, but it is the same 
for the other values). They have roughly the same minimum radial numbers n. (For m = 1, 
there are seen to be some with a ratio of 1.45.) However, these frequencies are part of an 
effective continuum of radial mode numbers which extends to higher (and possibly somehat 
lower) frequencies, as also discussed below. We do not know any obvious way in which such 
a ratio could be observed. 

Another possibility to explain the observed ratio is related to a pair of low-frequency 
fundamental [m = j = 0) p-modes found in RD3. According to Table 1 of this paper, the 
ratio of the eigenfrequencies of such modes with radial mode numbers n = 2 and = is 
1.53. However, the capture zone of these modes seems too small to produce the observed 
amount of luminosity modulation. 

From this point of view, a better candidate is the m = 3 and m = 2 pair of g-modes. As 
found in RDl, their eigenfrequencies |ct| = mfi(rj) for m > 1, and their radial and vertical 
mode splitting is much smaller. Of course, it is not obvious why the m = and m = 1 
g-modes, as well as the j = —1 p-mode (Table 7) and the n = 1 low-frequency p-mode 
(RD3), are not also observed. However, the (turbulent?) excitation of modes is not well 
understood. 



We should note that lArras. Blaes. and Turner! (120061 ) found no g-modes in their shearing 



box MHD simulations of a limited radial region of an accretion disk. There were some 
indications of the generation of p-modes within the MRI-induced turbulence, however. 

The full-disk p-modes studied in this paper share some characteristics with the other 
modes studied previously. For example, their eigenfrequencies show a dependence on the 
properties of the central object and are rather insensitive to changes in the luminosity of the 
disk. Like most of the other modes, the properties of the full-disk p-modes are influenced 
by effects of general relativity in a strong field regime, and like the g-modes and some of the 
c-modes, they include the hottest part of the disk. 

On the other hand, the full-disk p-modes have some distinct characteristics. In addition 
to the fact that they cover the whole disk, there seem to be no full-disk p-modes with low 
or moderate values of the radial mode number n. We also note that in the outer parts of 
the disk, their eigenf unctions mainly occupy only a small layer close to the surfaces of the 
accretion disk. 

The radial wavelength of a minimum n mode is given by equation (18-11) as ~ (lOar)^^. 



-37- 



The observable effects of such modes are hard to perceive because their integrated luminosity 
fluctuations will average to effectively zero over the disk. It would seem that the most likely 
to be observed would be those axisymmetric modes whose frequency \a\ lies closest to the 
peak value of K,{r) (as was the case for the traveling waves mentioned above). In this case, 
the wavelength would be largest in that region, which is also located near the hottest part 
of the accretion disk. 

We now note some importa nt effects of viscosity on these (adiabati c) oscillations. If one 



applies the analysis developed in lOrtega- Rodriguez and Wagoner! (|2000| ) to the full-disk m = 
p-modes, one finds that they are viscously unstable, with a growth rate 1/t > a^a"^ /Q{ro) 
obtained approximately. (A similar result holds for the other modes with small values of m.) 
Here a* is the ratio of the shear to the total pressure, with the shear viscosity rj ~ a^h'^Q. 
This instability thus exists in both radiation and gas pre ssure dominated regions of the disk , 
unlike the axisymmetric Lightman-Eardley instability ( Kato. Fukue. and Mineshigel 1998 ). 



Our approach generalizes the local analysis of acoustic waves, which are also found to be 



viscously unstable under certain conditions (jKato. Fukue. and Mineshigdll998l ). 



Although this result violates the assumption that the imaginary part of the eigenfre- 
quency is small compared to the real part, it does suggest that the accretion disk is viscously 
unstable. This may well be a significant source of turbulence. In addition to the fact that 
the viscous width of the modes is much greater than the frequency separation of the radial 
harmonics, this growth would lead to the production of an effective continuum of excitation. 
We hope to investigate viscous damping and/or instability within a broader framework in the 
future, including a more realistic model of the effective (magnetically-generated) turbulent 
viscosity. 

Various p-modes were studied here in rather general terms, employing a WKB approach 
to solve the eigenvalue problem for the separated differential equations. As we have explained 
in Section 4, the nonexistence of some of these modes can be understood in physical terms 
taking into account the necessary continuity of the solutions together with the fact that 
formal conditions must be met over rather large radial regions containing changing physical 
conditions (namely, the transition from gas pressure to radiation pressure domination in the 
inner disk). 

On the other hand, there exists the possibility that some other modes do exist but were 
not found due to the limitations of our method. Those limitations follow from the character 
of the WKB approach, only valid in principle for sufficiently small wavelengths. Applied to 
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our study, this means that WKB is only strictly applicable when ^ 1, although there 
are indications that it works for Icj^,! ~ 1 as well. Having said this, it must be stressed that 
the modes found in these and the previous papers of the series provide a basis for dealing 
with the general problem of perturbations in thin accretion disks, assuming that the classes 
we have considered constitute a complete set. 
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We are grateful to Lev Kapitanski for his help with proving the spectrum estimate fl4-5p . 



Appendix. Lower Bound for Vertical Eigenvalues 

We introduce a new unknown function in the vertical equation (12-101) as 

v{y) = (1 - y'Y^Xiy), Vyiy) = (1 - yT'^My) , 

with some real 6. As Vy{y) and its derivative are finite at = 1, for v{y) and its derivative 
at j/^ ^ 1 — we obtain 

v{y) = 0{{1- yy/') , v'{y) = O {{1 - y^"'') . {A.2) 

Equation (12-101) in terms of v{ii) reads: 

(1 - y^) v" + 2y{g- 6) v - ai(l - y^) v + a^il - y\\ + \{l-y^)v = Q, (A.3) 

where 

A = 2gul {^/ul - 1) , (A.4) 

and 

ai(5, g, ul) = S^- {2g - 1)6 - 2gul = a2{5, g) - 5 ~ 2gul 

{A.5) 

a^iS, g)^6'- 2{g - 1)5 = {5 - g + 1)^ - {g - 1)^ . 

To make the following formulas shorter, we introduce some notations. Given any real 6 
and any function f{y) on y E (—1, 1), we set a one-parameter family of norms 

'\l-y'yf{y)dy. (A6) 
1 
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Since, for an arbitrary e > 0, 



we have an important inequality 
provided that the last norm, and hence all of them, is finite (evidently, \\f\\l < ||/||^/ for 

e>e'). 



Choose now some real 7 satisfying 

7 + 5 > , {A.8) 

multiply equation (A. 3) by (1 — y'^)'^v{y), and integrate over (—1, 1). Because of the above 
condition and formulas (A. 2), all the integrals converge and all the double substitutions at 
y = ±1 vanish in the following repeated integration by parts [we drop the limits ±1 in all 
the integrals and use the identity y^ = 1 — {I — y^) wherever necessary]: 

y"(l-|/2)^+Wdy = -||t;'||^+, + 2(7+l) J y{l - y'yvv' dy ; 
2 / y{l-yyvv'dy = -\\v\\^^ + 2^ [ y\l-y^y~-Wdy = 



- 1 1^ + 27 1 _ 27 1 |2 = -(27 + 1) 1 |2 + 27 I \v\\^^_^ . 
With these results, the described integration of the equation (A. 3) yields 

\ I I 1 12 1 1 ' 1 12 I 1 1 I |2 1 1 I |2 \ I I 1 12 1 1 I |2 f A c\\ 

-^IfI^+i = IP II7+1 +Pi IPII7 -P2 |p||-y_i > Pi\\v\\^ - P2\\v\\^_i 

for any eigenfunction v corresponding to the eigenvalue A [they are related to the eigenfunc- 
tion Vy and eigenvalue ^ by equahties (A.l) and (A.4), respectively]. Here 

Pi =Pi{l,^,g,^l) = (27 + 1)(7+ ^ + S - 9) + ai{S,g,(jjl) , 

(^.10) 

P2 = ^2(7, 5, g) = 27(7 + l + S-g) + 02(5, g) . 

Prom the point of getting the lower bound for A the problem with the inequahty (A. 9) 
is that the norm on the left is smaller than both norms on the right. One way to resolve it 
seems to be the following. 
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Suppose that for some given g > 0, lj^ > there exist 7 and S satisfying inequahty 
(A. 8) such that 

pi(7, 5, g, ul) < 0, p2(7, ^, ^7) < , (All) 

which makes the coefficient in front of the 7-norm in (A. 9) negative, while the other one, 
in front of the biggest (7 — l)-norm, positive. Then we can use inequahty (A. 7) with 
f = V, 9 = 'J to write (recall that pi is negative). 



II 11^ ^ II l|2 , Pi II 



and then set 

Pi 



> 

2P2 



to eliminate the biggest norm in the inequality (A. 9). The latter thus reduces to 

X\\v\ 



|2 > i^|b||2 



4P2 

or simply to 

A>i%%4^ (A12) 

[recall that {pl/p2) < by the assumption (A. 9)]. 

This is the desired lower bound of the spectrum. If moreover 11 = 'n.{g^uj1) is the set 
of all admissible values {7,5} for given g and ujI, then (A. 12) provides exactly the estimate 
(14-51) with the best constant 

A{g,uJl) = -- sup ^ = - mf - — r rr > . (A.13) 

4 {^,5}en P2[l,o,g,) 4 {7,5}en \p2[l,o,g,)\ 

The only thing that remains to be demonstrated is that the set n(5f, ujI) is not empty, 
in other words, that the inequalities (A. 8) and (A. 11) are always compatible for at least one 
pair {7, 5}. The help with solving those three inequalities, two of which are quadratic in both 
7 and 5, comes from a geometrical approach which turns out effective if we reparametrize 
the inequalities using 

x = 7 + 5 + l-5f (^-14) 
instead of 5. It is easy to see that expressions (A. 10) for pi and p2 can be now rewritten as 

pi = (7 + l/2)2 + a;2-r?, ^3 = 7^ + a;^ - '"/j , (A15) 
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with 

Ti = rj{g,ujl) = ^{g-ll2f + 2gujl m = rjj{g) = - 1| . (A.IQ) 

Therefore inequahties (A. 11) specify two open circles, Cj and Cu, in the plane {7, x}, while 
condition (A. 8) corresponds to the open half-plane H specified by a; > 1 — g. Everything 
depends now on how g compares to unity. 

Indeed, if 5- > 1 (1 < T < 2), then 



n = y/{g - 1/2)2 + 2gul > m + 1/2 = g - 1 + 1/2 = g - 1/2, 

so that the circle Cu lies completely in the circle Cj. It also belongs to the half-plane H, 
hence, in fact, this circle is the intersection of all the three domains. Therefore we have the 
exact lower bound of the eigenvalues (A. 13) with 

n = nig) = Cnig) = {{l,S} ^^ + (^ + S+l-gr < {g-lY} = {{^,x} : + < rj^} . 

It is now straightforward to see that the minimum value in (A. 13) is achieved when x = 0, 
so that 

A(^,c^) = - inf 2 >0- (^-17) 

4 |7l<m rjj-Y 

However, finding explicit value of A requires solving analytically a cubic equation with the 
parameters r/ and ru in the coefficients, which makes it impractical. On the other hand, any 
particular values of parameters from Cu can be introduced into (A. 12) to give the estimate 
with a somewhat worse but explicit constant. In particular, inequality fl4-6l) is obtained 
when X = Q{5 = g — 1) and 7 = 0, which corresponds to the center of the circle Cu. 

If, however, = 1 (F = 2), then the circle Cu shrinks into a point, p2 is nonnegative 
everywhere, and no estimate is available. Moreover, for < 1 (F > 2), the circle Cu is 
outside if, inequalities (A. 8) and (A. 11) are incompatible, 11 is empty, and no lower bound 
for the vertical eigenvalue spectrum can be gotten in this way. Thus the question whether 
such a bound exists at all when g <1{T > 2) remains so far unanswered. 

One could note that A would be nonnegative (no p-modes would exist) if > and 
< in the inequality (A. 9). The above analysis shows, nevertheless, that such situation 
is impossible for any g > 
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Table 1. Eigenfrequency boundary — (T_ = max[fi;(r) + mfl{r)]; fl{ri) (last row) 



a 



m 0.5 a± 0.953 0.998 



0.02210 0.03312 0.06597 0.07492 

1 0.07720 0.1209 0.2962 0.4277 

2 0.1414 0.2244 0.5658 0.8457 
nin) 0.06804 0.1086 0.2784 0.4213 



Table 2. Radii that maximize K{r) + mr2(r) 



a 



m 0.5 a± 0.998 



8.0000 5.7628 2.9639 2.1929 

1 6.4000 4.6377 2.0073 1.2608 

2 6.1407 4.3000 1.9497 1.2424 



Table 3. Relevant radii in the Novikov & Thorne disk solution, for M/Mq = 10 and 

L/L^-M = 0.1 



a Ti Ri R2 



6 6.92 140 

0.5 4.233 4.81 100 

0.953 1.9099 2.06 53 

0.998 1.23697 1.26 34 



Table 4. Values of 7r(j + 1 + g/2)/^/2g for even modes with (7 = 3/2 and g = ?> (gas and 
radiation-pressure dominated regions, respectivelly), relevant for condition (15-291) . 



3 



-1 


1.36 


1.92 





3.17 


3.21 


1 


4.99 


4.49 


2 


6.80 


5.77 


3 


8.62 


7.05 



Table 5. Values of i^J for even modes with = 3/2 and (7 = 3 (gas- and radiation-pressure 
dominated regions, respectivelly), relevant for the condition 1 — ^/cj^ > 1. 



3 






-1 


3 


6 





16.33 


16.67 


1 


40.33 


32.67 


2 


75 


54 


3 


120.33 


80.67 



-47- 



Table 6. Values of li (shorthand for lij). 



a 



To 


3 







0.5 


a_L 




_i 


r 
-'1 




4 '^fi 1 n^ 


fi '^fi 1 n^ 






-'2 


1 An 1 


1 49 1 n^ 


9 97 1 n^ 








94 ^ 


'^n 9 








-'4 


fin fi 


85 8 


SI n 




n 


r 


9 5fi 1 n"^ 


9 51 in^ 


77 1 n'^ 

O. 1 1 xu 








72 S 


754 


1 22 1 n^ 








57 1 


7n n 


195 






-'4 


29.2 


41.6 


149 


1 nn 


_1 


r 

-'1 


7 '?9 1 n^ 

/ .OZ lU 


Q 7^ 1 n^ 


4s 1 n''' 






^2 


1 "^f; 1 n^ 


1 7*^ 1 n^ 


^71 1 n^ 






^3 


oo.z 


4n Q 


■iQ 4 






T 

-'4 


oo.U 


izi 


/1 1 7 







/l 


4.39 10^ 


5.85 10^ 


1.52 10^ 






/2 


7.9710^ 


1.0310^ 


2.60 10^ 






h 


84.2 


104 


207 






h 


45.6 


62.8 


201 




-1 


h 


2.21 10^^ 


2.38 10^^ 


2.82 10^^ 






h 


4.51 10^ 


4.85 10^ 


5.75 10^ 






h 


-2.04 


-3.19 


7.05 






h 


417 


494 


857 







h 


9.49 10^2 


1.0210^^ 


1.21 10^^ 






h 


1.9310^ 


2.08 10^ 


2.4710^ 






h 


305 


352 


498 






h 


188 


223 


390 
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Table 7. Lowest values of the eigenfrcqucncy |a|, the radial number n and the critical 
radius Tc for full-disk modes with ^ < and m — 0. 



a 



' o 


n 

J 




n 






20 


-1 


\a\ 


0.117 


0.176 


0.322 






n 


205 


439 


2.12 10^ 








Ri 


Ri 


Ri 







\a\ 


0.238 


0.343 


0.537 






n 


479 


960 


3.50 10^ 












Ri 


100 


-1 




0.117 


0.176 


0.322 






n 


3.1710^ 


9.65 10^ 


1.15 10^ 






Tc 


Ri 


Ri 


Ri 







cr 


0.238 


0.343 


0.537 






n 


7.93 10^ 


2.1910^ 


1.40 10^ 












Ri 


10^ 


-1 




0.117 


0.176 


0.322 






n 


9.5710^° 


2.35 10" 


9.3310" 








Ri 


Ri 


Ri 







|(7 


0.238 


0.343 


0.537 






n 


1.7110" 


3.81 10" 


1.1110^2 












Ri 




1 


\(t\ 


0.374 


0.538 


0.842 






n 


2.6910" 


5.9810" 


1.7410^2 
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Table 8. Lowest values of the eigenfrequency \a\, the radial number n and the critical 
radius Tc for full-disk modes; ^<0, m = 1, a < case. 



a 





3 




n 




a± 


20 


-1 


a 


0.171 


0.267 


0.578 






n 


176 


827 


7.02 10^ 








Ri 


Ri 


Ri 







\a\ 


0.306 


0.451 


0.814 






n 


707 


1.65 10^ 


8.30 10^ 














100 


-1 




0.171 


0.267 


0.578 






n 


6.5710^ 


2.1810^ 


3.70 10^ 








Ri 


Ri 


Ri 







cr 


0.306 


0.451 


0.814 






n 


1.3010^ 


3.79 10^ 


3.22 10^ 














10^ 


-1 


\a\ 


0.171 


0.267 


0.578 






n 


2.0610^^ 


5.4010^^ 


3.00 10^2 








Ri 


Ri 


Ri 









0.306 


0.451 


0.814 






n 


2.8310^^ 


6.6010^^ 


2.55 10^2 
















1 


(7 


0.442 


0.647 


1.12 






n 


3.7610^^ 


8.65 10^^ 


3.08 10^2 
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Fig. 1. — Universal curve for the vertical WKB equation from the analytical asymptotic 
(top) and numerical (bottom) solutions. 
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Fig. 2. — Dependence of {(3-s/g^)ri on the angular momentum of the black hole. 
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